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Abstract 

The Hubbard model describes interacting electrons on a lattice, a 
situation which occurs in many solid state materials and devices. The 
aim of the present paper is to briefly discuss this model and its appli- 
cations in the study of transport properties of quasi-one dimensional 
and quasi-two dimensional systems. Several examples are taken from 
materials and fields intensively studied in recent years: organic con- 
ductors and semiconducting heterostructures. 



1 Introduction 



All natural sciences are devoted to an easily formulated aim: trying to dis- 
cover,explain and predict how natural phenomena occur. Performing exper- 
iments on any kind of a specimen in any natural science means putting the 
system under study into interaction with some kind of external probe,and 
then measuring the system's response. If this interaction is weak,the sys- 
tem's response is linear and in solid state physics this is usually called the 
one electron picture. 

It has been shown in the last few decades that there exist materials which 
can not be described within the one-electron picture. Examples of such ma- 
terials are the organic conductors and the high-temperature superconduc- 
tors. Quasi-one- dimensional (QID) organic conductors were synthetized in 
the last century [1]. Their general chemical formula is {TMTSF)2X where 
(TMTSF)2 is a comphcated compound called bi-tetra-methyl-tio-seleno- 
fulvalene,and X is any anion which can be added to the central complex. 
Some examples of the anions which can be added are: ClOi,FSO^,ReO^..¥oY 
a recent review and some history of the field see [2] . 

Early experiments have already shown that the electrical conductivity 
of these new materials,which later became known as the Bechgaard salts, can 
not be described within the standard theory of metallic conductivity [3] . This 
discrepancy was explained by the fact that the electrons in Bechgaard salts 
are correlated. 

The simplest theoretical model applicable to the Bechgaard salts is the 
Hubbard model (HM) [4]. The present paper is devoted to a brief intro- 
duction to this model and its transport properties, and applications of these 
theoretical results to two kinds of systems, under the limitation of low di- 
mensionality. This limitation is motivated both mathematically and physi- 
cally. Mathematically,the Hubbard model has so far been solved only for one 
dimensional systems (ID) [5]. Physically,it has been shown in recent years 
that low dimensional materials have extraordinary properties. 

2 The Hubbard model 

Around the middle of the last century Mott developed a model of the metal 
to insulator transition (MIT), and showed what infiuence correlations can 
have on the results of the one electron picture [6] . The developement of the 
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HM was motivated by the need to develop a microscopic model of the Mott 
transition. 

The starting ingredient of the HM is a collection of atomic orbitals.The 
model is developed under the assumption of tight binding. This means that 
an electron's wave function is assumed to be centered on one atomic site, and 
that an electron can "hop" just one lattice spacing at a time. The rigor- 
ous formulation of the HM turned out to be too complex for any practical 
applicability, so Hubbard had to introduce considerable simplifications. Ba- 
sicallythe Hamiltonian is given as a sum of a "free" kinetic term Hq and the 
interaction term Hj. It can be shown that the simplest expression for the 
Hubbard Hamiltonian in one spatial dimension is: 

H=-t J2i4+l,aCi,a + ci^aCi+l,a) + UY^ 71^,^^^,; (l) 
i,rT i 

where t is the hopping integral, [/ the on-site repulsion, and all the other 
symbols have their usual second quantisation meanings. 

To study the interplay of interactions and disorder in electronic systems 
disorder can be introduced into the Hubbard model. The most common type 
of disorder is that of a random potential; for the Hubard model this trans- 
lates into random site energies. An equivalent formulation is that of random 
bonds,which practically means random hopping parameters. If we allow for a 
paralel magnetic field which couples to the electron spin,the Hamiltonian for 
2D calculations becomes 

B denotes the magnetic field which couples to the electron spin, and n is 
the chemical potential. The strength of the disorder is tuned by taking the 
hopping parameters tij from a probability distribution of the form P{tij) — 
for tij e [t — At/2,t + At/2] and zero otherwise. 

3 The transport properties 

Studies of the conductivity of the Bechgaard salts are of interest only under 
low temperature conditions,which means that electron-electron scattering is 
the dominant scattering mechanism. The electrons in these materials can be 
described as a Fermi (or Luttinger) liquid. Accordingly, it should be possible 
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to apply methods of the theory of fluids to the calculation of the transport 
properties of the Bechgaard salts. Various methods for the calculation of the 
conductivity,assuming that the Hamiltonian is known, have recently been re- 
viewed in [7]. In work on the Bechgaard salts,the "memory function" method 
was used [7]. 
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This method was developed in the seventies, and within it the conductivity 
can be calculated as follows: 



Xab(z) =« A;B »= -i / e^izt < [A(t),Bm > dt (3) 

Jo 

where A — B — [j, H] ,j denotes the current operator, and 

a{z) = t^[l-^ (4) 



4 Results 



4.1 The Bechgaard salts 

The electrical conductivity of the Bechgaard salts was calculated using Eqs. 

(1),(3) and (4). Full details of the calculation have been discussed at 
length in [8]. It was assumed there that the electrons in a Bechgaard salt 
can be approximated as a Fermi hquid (FL). Field theory shows that the 
FL model fails in ID systems (for example [10]). However, the Bechgaard 
salts arc experimentally known to be quasi one-dimensional, although recent 
experiments suggest that a 2D Fermi surface exists in some of them [11], 
which can be taken as a justification for the applicability of the FL model. 

The first step in the calculation was to determine the chemical potential of 
the FL on a lattice. The resulting expression has the form of a quotient, and 
fits to the known result of Lieb and Wu [5]. The next step was the deter- 
mination of the susceptibility Xt^^'i the final one was the calculation of the 
conductivity. Various mathematical approximations were made, and they are 
discussed in [8]. 

The following figure represents an example of the temperature dependence 
of the conductivity of the Bechgaard salts. The curve is normalized to unity 
at T = 116K and it is drawn for two different values of the band filling n. 
Experimentally, varying the band filling means doping the material with an 
electron donor or an acceptor.The calculation was performed without taking 
into account the influence of high external pressure, which is usually varied 
in experiments (for example [9]). 
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Figure 1: Normalzed conductivity of a Bechgaard salt for two values of the 
band filling 



Later work has shown that a suitable change of variables gives the pos- 
sibility of taking into account the influence of high external pressure on the 
electrical conductivity [12]. The real part of the susceptibility has the form 

where qi are numerical constants. The functions Ai contain all the param- 
eters of the problem which means also the value of the lattice constant, except 
the frequency. Starting from this result, one can calculate the immaginary 
part of the susceptibility, and after that the real part of the conductivity. 
The final result for the real part of the conductivity has the following form 

1 A- H{^f] 



As the equation of state of the Bechgaard salts is not known, the influence of 
high pressure on the conductivity is taken into account through the change 
of the lattice constant. 

Apart the electrical conductivity,experiments in the past ~ 10 years have 
given the possibility of measuring the thermal conductivity of the organic 
conductors. First experiments were made on the Bechgaard salts, but these 
were quickly followed by studies of quasi two dimensional organic conductors. 
Examples are [16], 17]. 

4.2 Two dimensional systems 

The basic numerical methos used in applying the Hubbard model to two 
dimensional systems is the so called Determinant Quantum Monte Carlo 
Method (DQMCM), developed at the University of Santa Barbara. A brief re- 
view of the method is presented in [13] and references given there. Details on 
the calculation of the partition funaction are also available there. Interesting 
results have been obtained in applying DQMCM to the 2D Hubbard model 
[14]. Plots of the conductivity as a function of the temperature T for various 
values of the disorder show an indication of a metal to insulator transition 
above a critical disorder of Ac = 2.7. If such a behaviour persists to very low 
temperatures, and in the thermodynamical Imit, it would describe a ground 
state metal to insulator transition driven by disorder. The temperature is ex- 
pressed in units of t. The parameter varied in laboratory experiments is 
the electron density, while in theoretical work on the Hubbard model one 
varies the disorder strength A. It can be shown that the two approaches are 
equivalent. 

It has been shown [15] that a Zeman magnetic field reduces the conduc- 
tivity of a conducting disordered electronic 2D system with a fixed disorder 
and a variable field, and it becomes temperature independent after a certain 
critical value of the field. At larger field strenghts there are indications of 
insulating behaviour. The occurence of a metal to insulator transition driven 
by a magnetic field is in agreement with real laboratory experiments such as 
[18]. In that particular experiment, the value of the resistivity at the transi- 
tion was found to be density dependent, which implies that it can be controlled 
by high external pressure. It is interesting to note that a recent study of a 
2D Hubbard model with site disorder gives physically similar results [13]. 



6 



5 Conclusions 



In this paper we have briefly discussed the Hubbard model and some of its 
apphcations to the study of transport properties of quasi-one-dimensional 
and quasi-two-dimensional systems. Due to space limitations, many inter- 
esting aspects have not been discussed. In the ID case,the Hubbard model 
was apllied to the Bechgaard salts, at first without but then with taking into 
account effects of external pressure. 

Concerning possible metal to insulator transitions in 2D electronic dis- 
order systems, it has been shown that interactions enhance the conductivity. 
At low temperature, a transition from a metallic state at weak disorder to an 
insulating state at strong disorder can occur. 

In both dimensionalities,it would be interesting infuture work to allow for 
longer range interactions. It seems that a necessary step in applications of 
the Hubbard model is to investigate in detail its thermodynamcal properties 
like the compressibility or susceptibility. 
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